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. ' Abstract 

The nonsymmetric gravitational theory predicts an acceleration law that 
! modifies the Newtonian law of attraction between particles. For weak fields a 

^vO ' fit to the flat rotation curves of galaxies is obtained in terms of the mass (mass- 

to-light ratio M/L) of galaxies. The fits assume that the galaxies are not 
dominated by exotic dark matter. The equations of motion for test particles 
reduce for weak gravitational fields to the GR equations of motion and the 
predictions for the solar system and the binary pulsar PSR 1913+16 agree 
with the observations. The gravitational lensing of clusters of galaxies can be 
explained without exotic dark matter. 
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99; 1 Introduction 

> ■ 

^ . A gravitational theory explanation of the acceleration of the expansion of the uni- 

verse [H 13 E] and the observed flat rotation curves of galaxies was proposed [I] , 
based on the nonsymmetric gravitational theory (NGT) jSl El EI- Since no dark 
matter has been detected so far, it seems imperative to seek a possible modified 
gravitational theory that could explain the now large amount of data on galaxy 
rotation curves. The same holds true for the need to explain the acceleration of the 
expansion of the universe without having to invoke a cosmo logical constant, because 
of the serious problems related to this constant [H|. 

In the following, we summarize the derivation of the motion of test particles in 
NGT. We consider the derivation of test particle motion from the NGT conservation 
laws. The motion of a particle in a static, spherically symmetric gravitational field is 
derived, yielding the modified Newtonian law of motion for weak gravitational fields. 
Two parameters y/M and r occur in the generalized Newtonian acceleration law. 
The parameter \/M is modelled for a bound system by a dependence on the mean 
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orbital radius of a test particle, and the range parameter ro is determined for galaxies 
and clusters of galaxies from the acceleration cHq where Hq is the measured Hubble 
constant. A fit to both low surface brightness and high surface brightness galaxies is 
achieved in terms of the total galaxy mass M (or Mj L) without exotic dark matter. 
A satisfactory fit is achieved to the rotational velocity data generic to the elliptical 
galaxy NGC 3379. Fits to the data of the two spheroidal dwarf galaxies Fornax and 
Draco and the globular cluster u> Cenauri are also obtained. The predicted light 
bending and lensing can lead to agreement with galaxy cluster lensing observations. 

The modelled values of the parameter y/Mo for the solar system and Earth, lead 
to agreement with solar system observations, terrestrial gravitational experiments 
and the binary pulsar PSR 1913+16 observations. 

2 The Field Equations 

The nonsymmetric fundamental tensor g^ u is defined by 00 Ej: 

9nu = 9(jiu) + #[H> (!) 
where ^ ^ 

9Qa>) = r (</,«/ + SV)> SV] = nig ijv - 9vp)- (2) 

The nonsymmetric connection Y\ v is decomposed as 



— r ( ^) + r^. (3) 

The contravariant tensor g^ u is defined in terms of the equation 

g^g™ = g^g™ = 5V (4) 

The NGT action is given by 

^ngt — S + Sm, (5) 

where 

S = ^ / d"x[g^R; v (W) - 2Av^ - ^V^ M ], (6) 
and Sm is the matter action satisfying the relation 

1 W (7) 



Here, we have chosen units c = 1, = \/—gg^ w , g = Det(g M! ,), A is the cosmological 
constant, /z is a mass associated with the skew field g\p, v \- Moreover, is the 
nonsymmetric energy-momentum tensor and R*(W) is the tensor 

r; u (w) = r^(w)-Iw„w v , (8) 
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where R flu (W) is the NGT contracted curvature tensor 

RAW) = Wis - + W?) - W^W^ + W^W* (9) 



defined in terms of the unconstrained nonsymmetric connection 

2 
3 

where 



w£ = rjL, - o s \w v , (io) 



w, = l(w^-w^). (11) 
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Eq. (fTU|) leads to the result 



T, = rf M] = 0. (12) 
The contracted tensor R^ U (W) can be written as 

R,AW)=R, U (T) + ^W M , (13) 
where ^ 

R/j,u(f) = r^ uj3 — - (r^^ + rf^)^) — rf^r^g + r^r^. (14) 

The gravitational constant G in the action S is defined in terms of the "bare" 
gravitational constant Go: 

G = G Z, (15) 

where Go — 6.673 x 10~ 8 g _1 cm 3 s~ 2 is Newton's constant and Z depends on the 
strength of the coupling of g^ u ] to matter. Thus, Z = 1 when is zero and NGT 
reduces to Einstein's GR. 

A variation of the action S n g^ yields the field equations in the presence of matter 
sources 

G%(W) + Ag^ + = SttGT^, (16) 
g [H „ = ~\^ a) W a , (17) 

+^g w < /3] + \(&®w& - z m wM = 0. (18) 

Here, we have G*^(W) = R*^(W) - lg^K*(W), where TZ*(W) = g^R*^ u (W), and 

= \tf(g\M + \g^g wp] g[ P a\ + g [(Tp] g^g pv ). (19) 

The vacuum field equations in the absence of matter sources are given by 

R; U (W) = Ag, u - (S, u - ^S), (20) 
3 



where S = g^ u S^ v . For the case of a static, spherically symmetric field with \i = 
and <?[oi] = (i=l ,2,3) , the gravitational field is described by the Wyman solution 
given in Appendix A 9j. In this case the vector = jU |Hl E] and the field 
equations (|2Uj) in empty space T^ v = become 

R^(T) = A 9fM „ - (S^ - \g^S). (21) 
The time independent components of these field equations are given in Appendix A. 



3 The Equations of Motion of a Particle 



In NGT there are two choices for the motion of a particle [101, for there exist two 
connections, the Christoffel connection given by 



\s [xp) (g M ,u + g{p V )^ - gw, P ) , (22) 



where 

s {va) g {m) = r» (23) 

and the nonsymmetric connection If we define the parallel transport of a vector 
V by 

D,V x = d,V x + T x piM V p , (24) 
then we obtain the equation of motion 

u v D v u x = ^r + rJXu" = 0, (25) 

where r is the proper time along the path followed by the particle and u x = dx x /dr 
is the 4-velocity of the particle. This defines the path equation which is not a path 
of extremal length. 

We can derive the second equation of motion by using the Lagrangian defined 

by 

L = gwuW. (26) 

We have 

We see that when the right-hand side vanishes, the Euler-Lagrangian equations are 
satisfied and we get the extremal geodesic equation 
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The conservation laws in NGT can be written |TT] : 



d p f { - \d x g, v T^ = 0, (29) 

where 

f^ = i(^ A T^ + ^T^). (30) 
Let us assume a monopole test particle 

f d 3 xT^ u 0, (31) 

d 3 x(x a -X a )T^ u = 0, (32) 

and similarly for higher moments. The integration is carried out over a hypersurface 
of constant t, following the procedure of Papapetrou ^2]; and X a is the position of 
the monopole. It can be shown that (J2*9~|) leads to 

at \ at J J I \w I at at J 



d x g, v J d 3 xT^ - dxg { »^^- J rf 3 xt 00 ) , (33) 
where = s {Xu) f £ If we assume that T [ ^ ] = 0, then we get 

+m( M— — = (34) 
dr\ dr J \fj,v J dr dr 

where 

m= d ^ I d 3 xt 00 . (35) 



dt 

This leads us to the geodesic equation (j2HJ). However, we cannot in general assume 
that TlH is zero, for it is associated with the intrinsic spin of matter or to the 
skew field itself, and we must include a coupling to the spin and skew field on 
the right-hand side of ()34jl . We shall find that in the weak field approximation the 
path equation (f23|) and the geodesic equation (J2%|) become approximately the same, 
including a necessary coupling to a skew symmetric source. 



4 Particle Motion in a Static Spherically Sym- 
metric Gravitational Field 

It was shown in ref. ^D] that for large values of r in the static, spherically symmetric 
solution of the NGT field equations, the geodesic equation and the path equation yield 
similar physical results. Since we shall be concerned with the dynamics of galaxies, 
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we shall treat the case of the geodesic equation coupled to a skew symmetric source. 
We have 

^ + {^} uV = sW/M< (36) 

where 

( e r,auX \ 
f[afj,} = Xdia^—^H^xjg^J . (37) 

Here, the skew tensor ^ vXri is the Levi-Civita tensor density and H^ vX ] is given by 

H[^u\\ = -(<9a#[H + d^g[ u x] + d u gy\^), (38) 

and A is a coupling constant with the dimension of length that couples the skew 
field to the test particle. 

In the static, spherically symmetric field of NGT, the tensor H^ u x] has only one 
non- vanishing component 

H^] = ^d r gi 04>] = /' sin 9 (39) 



and we get 

d ( if 



(40) 



/a 7 (r 4 + / 2 ). 

The equations of motion for a test particle are given by 

d 2 r a' fdr\ 2 r/d9\ 2 /sin 2 9\ (d<j)\ 2 7' / dt\ 2 
dr 2 2a\drJ a\drj \ a J\drJ 2a\drJ 

+ d ( w \(*\ = 0i (41) 



a 7 ( r 4 + p)J \dr 



dr 

£i + t(±) + l± u x jL ) = 0, (42) 

dr 2 •y\drJ\dTj 7 dr \, AWr* + f 2 V 



7 \dTj \dTj 7^ Vy , a 7 (r 4 + / 2 ) 

d 2 9 2(d0\(dr\ n Jd(j)\ 2 
— 7 + - — — - sin 9 cos 9 — 



^ + ;UJU- si "UJ =°- (43) 



(ir 2 r \dr J \dr , 

The orbit of the test particle can be shown to lie in a plane and by an appropriate 
choice of axes, we can make 9 = ir/2. Integrating Eq.(jHJ) gives 

r — = J, (45) 
dr 



6 



where J is the conserved orbital angular momentum. Integration of Eq. (|42)l gives 



dt 
fa 



7 



At/' 



V / «7(r 4 + / 2 



(46) 



where i£ is the constant energy per unit mass. 

By substituting (f4l)j) into (pfTj) and using ([45)1. we obtain 



<i 2 r a' / dr\ 2 
fa 2 + 2^\fa) 



J 2 7 
— ? + 



A7/' 



1 2 



«7(r 4 + f 2 ) 



+ E 



1 d 

7 (ir 



A7/' 



a7 ( r 4 + J2) 



(47) 

Let us now make the approximations that Xf'/r 2 C 1, //r 2 C 1 and the slow 
motion approximation dr/dt <C 1. Then, for material particles we set E = 1 and 
(|4*7|) becomes 

d 2 r J 2 N + GM^ x d fr 



X- 



dt 2 r 3 r 2 c?r Vi 
where is the Newtonian orbital angular momentum. 
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5 Linear Weak Field Approximation 

We expand about a Ricci-flat GR background 

9 l » = 9<j2) + h l » + 0{h 2 ), (49) 

where <?£^ is the GR background metric. The skew field /i[ M „] obeys the linearized 
equation of motion in the GR background geometry 

V CT i^ CT + 4h [a0i B Pltav + n 2 h [H = 0, (50) 

where 

F[fMu\] = dxh^v] + d^h^x] + d u h[ XfJ ] (51) 

and V A and B Pimv denote the background GR covariant derivative and curvature 
tensor, respectively. Moreover, 

W„ = -2V A % A] . (52) 

For flat Minkowski spacetime, Eq. (j50|) reduces exactly to the massive Kalb-Ramond- 
Proca equation which is free of ghost pole instabilities with a positive Hamil- 
tonian bounded from below. 

For the static, spherically symmetric spacetime, the linearized equation of motion 
f!50|) on a Schwarzschild background takes the form 




where we assume that 7(7-) ~ l/a(r) ~ 1 — IGMjr and j3(r) = r 2 . The solution to 
this equation in leading order is jSl 13 CS] 



sG 2 M 2 exp(-fxr) ( GM 



/(r)= 3 (H ^ (1+/ir+ 



2 + fir exp(2 fir) Ei(l, 2 fir) (fir — 1) 



, (54) 



where Ei is the exponential integral function 



E l (n,x)=fdt eM tn Xt) . (55) 

The constant sG 2 M 2 /3 is fixed from the exact Wyman solution with A = by 
taking the limit fi — > 0. 
For large r we obtain 

l g G 2 M 2 exp(-/xr)(l+/ir) 
= 3 (^ ' (56) 

and for fiGM <C 1, we have (fir)^ ~ 1, giving 

/( r ) = i s G 2 M 2 exp(- / ur)(l +/ir). (57) 

This is a solution to the equation 

f"(r) - -f'(r) - fi 2 f{r) = 0. (58) 
r 

If we consider the expansion about a GR Schwarzschild background, then (J58j) is 
valid for fi 2 > 8GM/r 3 . 

To the order of weak field approximation, we obtain from Eq. (j48J) : 

d 2 r J 2 N _ GM + vay£ri 



dt 2 r 3 r 2 r 2 

where the constant a is given by 

\sG 2 M 2 fi 2 



(60) 



Here, A and s denote the coupling strengths of the test particle and the skew g[ 23 ] 
field, respectively. In Eq. ()59|). we have required that the additional NGT accelera- 
tion on the right-hand side is a repulsive force. This is in keeping with the weak field 
approximation result obtained in jSJE], which corresponds to a skew force produced 
by a massive axial vector spin 1 + boson exchange. 
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6 Orbital Equation of Motion 

Let us write the line element as 

ds 2 = -fdt 2 - adr 2 - r 2 (d6 2 + sin 2 Odcj) 2 ). (61) 

We set 9 = n/2 and divide the resulting expression by dr 2 and use Eqs. fl45|) and 
flU to obtain 



- 



2 T2 



( dr\ 2 J 



\drj ar 2 cry 



A7f 



V^tF+7 2 ) 



E 



2 



E 



a 



(62) 



We have ds 2 = Edr 2 , so that dsjdr is a constant. For material particles E > and 
for massless photons i£ = 0. 

Let us set u = 1/r and by using (plo - )). we have dr/dr = —Jdu/dcp. Substituting 
this into (p32*|) . we obtain 



Y= — 



^(r 4 + f 2 ) 



E 



ar 2 a J 2 



(63) 



By substituting dr/dcj) = —(l/u 2 )du/d(j) into (jBHJ) and using the approximation 7 ~ 
I/a ~ 1 — 2GM/r and //r 2 <C 1 and Xf'/r 2 <C 1, we get after some manipulation 

d 2 u EGM EXsG 2 M 2 ( 1 \/ 1 \ om 2 

where r = 

For material test particles = 1 and we obtain 

d 2 u GM „ %r 2 K ( \\( 1\ 

_ +1I= _ +3 GM^--exp(--)(l + -), (65) 

where .fT = XsG 2 M 2 /3?"q. On the other hand, for massless photons ds 2 = and 
E = and (JMj) gives 

d 2 v 

^-+u = 3GMu 2 . (66) 
a</> 2 

7 Galaxy Rotational Velocity Curves 

A possible explanation of the galactic rotational velocity curves problem has been 
obtained in NGT ^BJ- From the radial acceleration derived from (fKHJ) experienced 
by a test particle in a static, spherically symmetric gravitational field due to a point 
source, we obtain 

a ( r ) = — + a 2 1 + — ' ( 6T ) 
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Moreover, G = Goo is defined to be the gravitational constant at infinity 



Goo = G (l + ^), (68) 

where Go is Newton's "bare" gravitational constant. This conforms with our defi- 
nition of G in Eq. (jl5|) . which requires that G be renormalized in order to guarantee 
that (jBTf) reduces to the Newtonian acceleration 

^Newton — \yd) 

at small distances r<r . The constant a is given by 

° = ^rf (?0) 

The integration constant s in (fTUj) . occurring in the static, spherically symmetric 
solution (see, Appendix A), is dimensionless and can be modelled as 

s = gM\ (71) 

where M is the total mass of the particle source, g is a coupling constant and a is a 
dimensionless constant. We choose a = —3/2 and XgGl/3c 2 rl = G \/M where M 
is a parameter. The choice of a = —3/2 yields for a galaxy dynamics the Tully-Fisher 
law [HI. 

We obtain the acceleration on a point particle 

^ = -^ + «„^^M( 1 + L). (72) 

By using (jo^|) . we can write the NGT acceleration in the form 



fl(r) = -— 1 1 + Vm 



T 

1 - exp(-r/r ) ( 1 H 

r 



(73) 



We conclude that the gravitational constant can be different at small and large 
distance scales depending on the size of the parameter a/M . We have two parame- 
ters: the parameter ^M and the distance range r . We assume that G^ scales for 
constant M with increasing strength as y/Mo, while for fixed values of M we have 
Coo - * Gq as the total mass of the source M —>■ oo. Let us model the parameter 
a = \/Mq as a function of the mean orbital radius of a test particle 

a= jM~ Q = k(r olh ) n , (74) 

where k and n are constants. We shall choose the exponent to be n = 3/2. We 
shall apply (jTSj) to explain the flatness of rotation curves of galaxies, as well as the 
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approximate Tully-Fisher law [19], GqM ~ t> 4 , where v is the rotational velocity of 
a galaxy and M is the galaxy mass 



M M, ■ M m + M DB + M f . (75) 

Here, M*, Mm, M^b and Mf denote the visible mass, the mass of neutral hydrogen, 
possible dark baryon mass and gas, and the mass from the skew field energy density, 
respectively. In |4 3 , we obtained from the modified Friedmann equations 

Q = Q b + Q m + Q f , (76) 

where Q b , Qm and Qf denote the fractional values of baryons, field matter 
density and smooth j? , field density background, respectively, obtained from the 
expansion 

#[H = g[„u] + <%H> ( 77 ) 

where the smooth background gf^< and the fluctuations Sg^ describe the "dark 
energy" and "dark matter" components, respectively. The NGT cosmological field 
equations must be solved to yield the observational results Q b ~ 0.05, Q m ~ 0.3, 
and Q ~ 0.7 in order to be consistent with WMAP and supernovae data [H El Ej. 
The quantities Vt m and Qf in NGT replace the dark matter and dark energy. The 
mass Mf obtained from the skew field density p m is expected to contribute to the 
total mass M of galaxies. 

The rotational velocity of a star v is given by 



lG n M f /M r . i \ ( rMlVa 



Let us postulate that the parameters M and r give the magnitude of the con- 
stant acceleration 

G M 

a = —2-. (79) 
r o 

We assume that for galaxies and clusters of galaxies this acceleration is determined 
by 

a = cH . (80) 

Here, Hq is the current measured Hubble constant Ho = 100 hkm s -1 Mpc -1 where 
h = (0.71 ± 0.07) PI]. This gives 

a = 6.90 x 10~ 8 cms" 2 . (81) 

A good fit to low surface brightness and high surface brightness galaxy data is 
obtained with the parameters 

M = 9.60 x 1O U M , r = 13.92 kpc = 4.30 x 10 22 cm (82) 
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and M (or the mass-to-light ratio M/L). Substituting M = 9.60 x 10 n M Q into (f79"jl 
yields the tq in (j82j) . Thus, we fit the galaxy rotation curve data with one parameter 
M and the total galaxy mass M. Since we are using an equation of motion for point 
particle sources, we are unable to fit the cores of galaxies. A possible model for the 
galaxy cores is to assume for a radius r <C r c , where r c is the core radius, a rotation 
curve of an isothermal sphere in the ideal case where we can consider a massless 
disk embedded in it. Then, for r <C r c : 



W) ~ — g — J r, (83) 

where p c is the core density. For r ^> r c , the rotational velocity curve will be 
described by the NGT model ()78j) . Further investigation of this issue will require 
solving the field equations of NGT for a core mass density profile and will be the 
subject of future research. 

We can now fix the constant k in Eq. (J74|) by using the relation 



* = T^W* (84) 

V orb/9 



where a g and (r or b) 9 denote the values of these quantities for galaxies, respectively. 
For the mean value (r or b) 9 we choose 

(r orb ) 9 = 200kpc (85) 

and using the value of M in (|H2*j) we obtain 

k = 2.02 x 10~ 13 g 1/2 cm" 1 . (86) 

The fits to the galaxy rotation curves v in km/s versus the galaxy radius r in 
kps are shown in Fig. 1. The data are obtained from ref. [THj . 

In Fig. 2, fits to two dwarf galaxies (dSph) are shown. We assume that the 
relation between the velocity dispersion a and the rotational velocity v takes the 
simple form in e.g. an isothermal sphere model for which v ~ V2a. The error bars 
on the data [20] for the velocity dispersions are large, and in the case of Draco, due 
to the small radial range 0.1 kpc < r < 0.6 kpc, the Newtonian curve for 



*7) 



cannot be distinguished within the errors from the NGT prediction. However, it 
is noted that the NGT prediction for v appears to flatten out as r increases. For 
Draco M/L = 28.93 + 50.30(9. 58)(M /L ), whereas for Fornax M/L = 1.79 + 
O.72(O.4O)(M /L Q ). There is also an expected large error in the distance estimates 
to the dSph. Another serious potential source of error is that it is assumed that 
dSph galaxies are in dynamical equilibrium. The two studied here are members of 
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the Local Group and exist in the gravitational field of a larger galaxy, the Milky 
Way. Thus, the tidal interactions with the larger galaxy are expected to affect the 
dynamics of dSph galaxies and the interpretations of velocity dispersions [21]. These 
issues and others for dSph galaxies are critically considered in the context of dark 
matter models by Kormendy and Freeman [2*2] . 

We have also included a fit to the elliptical galaxy NGC 3379. The elliptical 
galaxy NGC 3379 has been the source of controversy recently [2S]. The velocities of 
elliptical galaxies are randomly distributed in the galaxy. However, the gravitational 
potential that would be experienced by a test particle star or planetary nebula in 
circular rotation about the center of the galaxy can be extracted from the line- 
of-sight velocity dispersion profiles. The data for R/R e s > 0.5 refer to planetary 
nebula. 

We use the mean values of the extracted rotational velocities for NGC 3379, ob- 
tained by Romanowsky et al. and find that the NGT predicted rotational veloci- 
ties agree well with their data. According to Romanowsky et al. there appears to be 
a dearth of dark matter in the elliptical galaxy which needs to be explained by dark 
matter models and N-body cosmological simulations. For Milgrom's MOND [23JI2E] 
it is argued by Sanders and Milgrom [26J that NGC 3379 is marginally within the 
MOND regime with an acceleration a ~ (a )Miigrom ~ 1-2 x lCT 8 cms~ 2 , so MOND 
should not apply to the elliptical galaxies. As we see from the fit to the data, the 
NGT results agree well with the data. Romanowsky et al. also give data for the 
two elliptical galaxies NGC 821 and NGC 4494, but the intrinsic circular velocities 
associated with the line-of-sight velocity dispersion profile data are not given by the 
authors, although the trends of the data are similar to NGC 3379. 

A fit to the data for the globular cluster to Centauri is shown in Fig. 3. The 
data is from McLaughlin and Meylan [27] . We use the velocity dispersion data and 
assume that the data is close to the rotational velocity curves associated with the 
velocity dispersion a p i.e. the isothermal sphere model relation v ~ V2a p holds. The 
fit to the data reveals that the predicted rotational velocity cannot be distinguished 
from the Newtonian-Kepler circular velocity curve within the orbital radius of the 
data. The authors conclude that there appears to be no room for Milgrom's MOND 
or dark matter, whereas the NGT results agree well with the data. 

In Fig. 4, we display a 3-dimensional plot of v versus the range of distance 
0.1 kpc < r < lOkpc and the range of total galaxy mass M used in the fitting 
of rotational velocity data. The red surface shows the Newtonian values of the 
rotational velocity v , while the black surface displays the NGT prediction for v. 

Table 1, displays the values of the total mass M used to fit the galaxies and the 
mass-to-light ratios M/L estimated from the data in references given in [TH] . 

There are now about 100 galaxies with available rotational velocity data [T%] . 
However, some of these galaxies are not well described by spherically symmetric 
halos or have some other disturbing physical feature, so we are unable to obtain fits 
to these data. 
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In Milgrom's phenomenological model j2lj ISHj we have 

V 4 = Go-^(«o)Milgrom, (88) 

where 

(ao)Miigrom = 1.2 x l(T 8 cms~ 2 . (89) 

We see that (|88j) predicts that the rotational velocity is constant out to an infinite 
range and the rotational velocity does not depend on a distance scale, but on the 
magnitude of the acceleration (a )Miigrom- in contrast, the NGT acceleration formula 
does depend on the radius r and the distance scale tq which for galaxies is fixed by 
the formula (fKUj) . 



8 Local and Solar System Observations 

We obtain from Eq. (|f)5jl the orbit equation 



d 2 u GM K . . . 



r 



3GM , , . 

+ ^^u\ (90) 



c 2 



where now K = XsG 2 M 2 /3c 4 rQ. Using the large r weak field approximation, and 
the expansion 

exp(-r/r ) = 1 - ^ + ^) + ... (91) 
we obtain the orbit equation for r < tq: 

d 2 u GM K GM 2 

^ + M = ^-J| + 3 — M ' (92) 



We can write this as 



where 



d 2 u „ GM 



+ u = N + 3—u 2 , (93) 



dcj) 2 c 



N=™-^ (94) 

We can solve Eq. (j93|) by perturbation theory and find for the perihelion advance 
of a planetary orbit 

Aw = ^(GM 9 -c%), (95) 

where Jn = (GMqp/c 2 ) 1 ^ 2 , p = a(l — e 2 ) and a and e denote the semimajor axis 
and the eccentricity of the planetary orbit, respectively. 

We choose for the mean orbital radius of the solar planetary system 

(r O rb) = 1 a.u. = 1.49 x 10 13 cm, (96) 
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which yields using (|71j) and k = 2 x 10 13 g 1//2 cm 1 : 



= ( v / Mo) = 1.2xlO 7 g 1 / 2 (97) 



and 



This result gives 



^) _ = 2.6 x 10- 10 . (98) 
AG 




10- 10 . (99) 

and G = Gq within the experimental errors for the measurement of Newton's con- 
stant G . 

We choose for the solar system K Q <C 1.5 km and use G = G^ = G to obtain 
from (|95|) a perihelion advance of Mercury in agreement with GR. 

For terrestrial experiments and orbits of satellites, we choose the mean orbital 
radius equal to the earth-moon distance, (r or b)e = 3.57 x 10 10 cm, which gives 

1.7 x HT 11 . (100) 

This also yields G = Gq within the experimental errors. 

For the binary pulsar PSR 1913+16 the formula (|95p can be adapted to the 
periastron shift of a binary system. Combining this with the NGT gravitational wave 
radiation formula, which will approximate closely the GR formula, we can obtain 
agreement with the observations for the binary pulsar. We choose the mean orbital 
radius equal to the projected semi-major axis of the binary, (r orb ) N = 7 x 10 10 cm, 

giving (\JMq/M)n = 5 x 10~ 14 . Thus, G = Go within the experimental errors 
and agreement with the binary pulsar data for the periastron shift is obtained for 
K N < 4.2 km. 

For a massless photon E — and we have 



d 2 u GM 



+ u = 3 u 2 . (101) 



For the solar system using ()99)1 this gives the light deflection: 



in agreement with GR. 



9 Galaxy Clusters and Lensing 

We can assess the existence of dark matter of galaxies and clusters of galaxies in 
two independent ways: from the dynamical behavior of test particles through the 
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study of extended rotation curves of galaxies, and from the deflection and focusing 
of electromagnetic radiation, e.g., gravitational lensing of clusters of galaxies. The 
light deflection by gravitational fields is a relativistic effect, so the second approach 
provides a way to test the relativistic effects of gravitation at the extra-galactic level. 
It has been shown that for conformal tensor-scalar gravity theories the bending of 
light is either the same or can even be weaker than predicted by GR j2Hl EH| • To rem- 
edy this problem, Bekenstein j2Sj has recently formulated a relativistic description 
of Milgrom's MOND model, including a time-like vector field as well as two scalar 
fields within a GR metric scenario. However, the time-like vector field violates local 
Lorentz invariance and requires preferred frames of reference. 

The bending angle of a light ray as it passes near a massive system along an 
approximately straight path is given to lowest order in v 2 /c 2 by 



J \a L \dz, (103) 



where _L denotes the perpendicular component to the ray's direction, and dz is the 
element of length along the ray and a denotes the acceleration. The best evidence 
in favor of dark matter lensing is the observed luminous arcs seen in the central 
regions of rich galaxy clusters |22J. The cluster velocity dispersion predicted by the 
observed arcs is consistent within errors with the observed velocity dispersion of 
the cluster galaxies. This points to a consistency between the virial mass and the 
lensing mass, which favors the existence of dark matter. 
From (|1U1|) . we obtain the light deflection 

where 



M = M ( 1 + V^)' (105) 

We obtain from Eq.([7I]) and k = 2 x 10~ 13 for a mean orbital cluster radius (r or b) c i ~ 

2 Mpc: 

a cl = ( V / Mo)ci = 3.1 x 10 24 g 1/2 . (106) 

For a cluster of mass M c \ ~ 10 13 M Q , we obtain 



IMo\ 
M ) 



~ 22. (107) 

cl 



We see that M ~ 22M and we can explain the increase in the light bending without 
exotic dark matter. 

From the formula Eq. ()72|) for r ^> tq ~ 14 kpc we get 

<*•) = (108) 
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We expect to obtain from this result a satisfactory description of lensing phenomena 
using Eq.flmSJ). 

The scaling by the parameter «j = (^M )i caused by the varying strength of 
the coupling of the skew field g[23] (r) = f(r) sin9 to matter due to the renormalized 
gravitational constant is seen to play an important role in describing consistently the 
solar system and the galaxy and cluster dynamics, without the postulate of exotic 
dark matter. 



10 Conclusions 

There is a large enough sample of galaxy data which fits our predicted NGT ac- 
celeration law to warrant taking seriously the proposal that NGT can explain the 
flat rotational velocity curves of galaxies without exotic dark matter. We do predict 
that there will be galaxy matter additional to that due to visible stars and baryons, 
associated with the energy density p m residing in the skew field g\pv]- It is interesting 
to note that we can fit the rotational velocity data of galaxies in the distance range 
0.02 kpc < r < 70kpc and in the mass range 1O 5 M < M < 1O U M . without 
exotic dark matter halos. We are required to investigate further the behavior of the 
NGT predictions for distances approaching the cores of galaxies, using a disk profile 
density. The lensing of clusters can also be explained by the theory without exotic 
dark matter in cluster halos. 

We are able to obtain agreement with the observations in the solar system and 
terrestrial gravitational experiments for suitable values of the parameter Mq. This 
required that we scale G and the parameter \/M as functions of the mean orbital 
radius (r orb ) of bound systems with the behavior «j = (v / Mo)i = (?"orb)^ 2 - 

A numerical solution of the NGT field equations for cosmology must be imple- 
mented to see whether the theory can account for the large scale structure of the 
universe and account for galaxy formation and big bang nucleosynthesis, without 
requiring the existence of exotic, undetected dark matter and a positive cosmological 
constant to describe dark energy. 



11 Appendix A: The Static Spherically Symmet- 
ric Solution and the Vacuum Field Equations 

In the case of a spherically symmetric static field [5], the canonical form of g^ u in 
NGT is given by 





—a 








w \ 







-P 


fsinO 










—fsm8 


-{3sm 2 6 





V 


— w 








7/ 



(109) 
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where a, [3, 7 and w are functions of r. The tensor g^ v has the components: 



( -V— 

1 w-'—a'y 



V 







w A — cry 



/3 2 +/ 2 
/CSC6> 
"/3 2 + f 2 







/csce 

/3 2 +/ 2 
/3CSC 2 6> 

/3 2 +/ 2 





W^—OTf 





110) 



_2 / 

2 — rv^ / 



For the theory in which there is no NGT magnetic monopole charge, we have w = 
and only the g[ 2 3] component of g^ u ] survives. 

The time independent field equations, Eq.(J2H), in empty space are given by 



I-, 1/N2 , , D21 , a' A' , i_fct_ 



-A — -[{Af + AB ] + — — + 



7 



2 7 



1 



—Act + -/! 



4a 



«/ 2 



27 \2a 2 7 



4" /? 2 + f 



2 • 



111^ 



(112) 



1 



^ 22 (r) = ^ S3 (r)co S e^ = - + -^ 

B(fA' + 2(3B) 



1 \(2fB-pA!\' 2fB-(3A' 



4a 



(a'j+j'a) (113) 



+ 



4a/5 



-A \x 



1 o f 2 



4^ /3 2 + / 2 ' 



^oo(r) 



2a 



4^ /3 2 + / 2 ' 



i2 [23] (r) = sin# 



-A 

2a\2a 2 7 2 

i2 [10 ](r) = o. 
R m (r) = o 



5a 



Af--li 2 f 1 + 



/? 2 



/? 2 + f 



V a 7 
sin#. 



(114) 

(115) 

(116) 
(117) 

(118) 
(119) 



Here, prime denotes differentiation with respect to r, and we have used the notation 

A = \n(f3 2 + f 2 ), (120) 
f? ~ Pf 



B 



'121 



P 2 + f 2 ' 

Let us assume the long-range approximation for which the fx 2 contributions in 
the vacuum field equations (J21j) can be neglected and that /i -1 > 2M. We then 
obtain the static, spherically symmetric Wyman solution for A = jHj: 



a 



7 = exp(z/j, 

(f 2 + P 2 )(Y) 2 

4M 2 7 
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(122) 
(123) 



M ? (l + i8 )exp(-,) sech 1 + ftv + 
A + i 2 
where Mi, A, i? and s are integration constants and z/ is an arbitrary function of r. 

We shall restrict our attention to the asymptotic condition guv) — > i]^ u where r]^ 
is the Minkowski flat metric tensor rj^ = diag(l, — 1, — 1, —1). This requires that 

sinh 2 B = -l, (125) 

so that we obtain 

M, 2 (l + is) exp(— v) , 9r 1 . . 
f + i P = A + i [ ~2 ( + ^ ] ' ( } 

This form of / + i/3 does not place any a priori boundary condition on g^u]- 
We now obtain the form of the static Wyman solution: 

7 = exp(z/), (127) 

M 2 (z/) 2 exp(-z/)(l + s 2 ) r . . „ X1 2 . . 

a = 1V ; ^ V — -[cosh(az/) - cos(H] , (128) 

2M 2 

f = 1 + ^ 2 exp(-z/)[(l - As)sinh(az/)sin(6z/) 



+s(l — cosh(az/)cos(5z/)][cosh(az/) — cos(bu)] , (129) 



where 

fl= (— 2 
and z/ is implicitly determined by the equation 

r 2 (l + A 2 ) 

exp(z/)(cosh(az/) — cos(5z/)) 2 — — (1 — As) [cosh (az/) cos (5z/) — 1] 

+ (A + s)sinh(az/)sin(&z/). (131) 
At large distances r — > oo, we have 

2Mi 1 / 1 \ , , 

° = 1+ (T+^F + b> (132) 

7=l-^ + 0f'], (133, 



f = -Ar 2 + l s Mf + o(^j. (134) 

If we assume that A = 0, then M = Mi where M denotes the mass of the particle 
source. 
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Table 1. Values of the total galaxy mass M used to fit rotational velocity data. 
Also shown are the mass-to-light-ratios Mj L with L obtained from ref . ^H] • 
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Fig. 1 - Fits to low-surface-brightness and high-surface-brightness spiral galaxy 
data. The black curve is the rotational velocity v versus r obtained from the 
modified Newtonian acceleration, while the green curve shows the Newtonian 
rotational velocity v versus r. The data are shown as red circles. 
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Fig. 2 - Fits to the data of two dwarf galaxies Fornax and Draco. The simple 
relation V ~ \/2a is assumed between the velocity dispersion a and the rotational 
velocity v. The black curve is the rotational velocity v versus r obtained from the 
modified Newtonian acceleration, while the green curve shows the Newtonian 
rotational velocity v versus r. The data are shown as red circles and the errors 
(not shown) are large and for Draco the Newtonian fit cannot be distinguished 
from the NGT fit to the data within the errors. 
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Fig. 3 - Fits to the elliptical galaxy NGC 3379 and the globular cluster u 
Centauri. The black curve is the rotational velocity v versus r obtained from the 
modified NGT acceleration, while the green curve is the Newtonian-Kepler velocity 
curve. For uj Centauri the black and green curves cannot be distinguished from one 

another. 




Fig. 4 - 3-dimensional plot of v versus the range of distance 0.1 kpc < r < lOkpc 
and the range of galaxy mass 5 x 10 6 M® < M < 2.5 x 10 11 M . The red surface 
shows the Newtonian values of the rotational velocity v, while the black surface 

displays the NGT prediction for v. 
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